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ABSTRACT 

The importance of shocks in nova explosions has been highlighted by Fermi’s discov¬ 
ery of y-ray producing novae. Over three years of multi-band VLA radio observations of the 
2010 nova V1723 Aql show that shocks between fast and slow flows within the ejecta led 
to the acceleration of particles and the production of synchrotron radiation. Soon after the 
start of the eruption, shocks in the ejecta produced an unexpected radio flare, resulting in a 
multi-peaked radio light curve. The emission eventually became consistent with an expanding 
thermal remnant with mass 2 x 10 and temperature 10"^ K. However, during the first two 
months, the >10^ K brightness temperature at low frequencies was too high to be due to ther¬ 
mal emission from the small amount of X-ray producing shock-heated gas. Radio imaging 
showed structures with velocities of 400 km s“* (d/6kpc) in the plane of the sky, perpendic¬ 
ular to a more elongated 1500 km s * (d/6kpc) flow. The morpho-kinematic structure of the 
ejecta from V1723 Aql appears similar to nova V959 Mon, where collisions between a slow 
torus and a faster flow collimated the fast flow and gave rise to y-ray producing shocks. Op¬ 
tical spectroscopy and X-ray observations of V1723 Aql during the radio flare are consistent 
with this picture. Our observations support the idea that shocks in novae occur when a fast 
flow collides with a slow collimating torus. Such shocks could be responsible for hard X-ray 
emission, y-ray production, and double-peaked radio light curves from some classical novae. 

Key words: novae, cataclysmic variables - binaries: general - stars: variables: general - 
stars: winds, outflows - white dwarfs - radio continuum: stars 


1 INTRODUCTION 

The most common type of stellar explosion, a nova occurs 
when accretion of material onto a white dwarf (WD) from its bi¬ 
nary companion ignites a thermonuclear runaway, expelling matter 
and causing the system to brighten. Because novae occur frequently 
- an estimated 35 yr“' in our galaxy alone (Shatter 1997) - they 
can be used to probe stellar explosions, astrophysical outflows, nu¬ 
clear burning, and accretion. However, to find an accurate model 
of these eruptions, one must account for the presence of shocks, 
bipolar outflows, rings, clumping, and other asymmetries that have 
been shown to be present in different nova systems (e.g Slavin et 


al. 1995; Shara et al. 1997; Sokoloski et al. 2008; Chomiuk et al. 
2014). 

Radio observations of novae furnish a useful tool for record¬ 
ing the evolution of a nova’s expanding shell and understanding 
the properties of the ejecta. In the radio, the ejecta remain op¬ 
tically thick to thermal emission at significantly lower densities 
than at shorter wavelengths. The ejecta can therefore typically be 
detected for months to years in the radio. In classical novae, ra¬ 
dio emission is generally observed to be dominated by thermal 
bremsstrahlung radiation, interpreted as emanating from an ex¬ 
panding shell (Hjellming 1979; Seaquist & Palimaka 1977; Bode 
& Evans 2008). 
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While non-thermal emission has long been evident in nova 
systems with red-giant companions, such as RS Oph (Taylor et al. 
1989), in classical novae with main sequence companions, the dom¬ 
inant features in historical radio light curves and spectra are gener¬ 
ally consistent with an expanding thermal shell producing free-free 
emission (Hjellming 1979; Seaquist et al. 1980; Hjellming 1996; 
Bode & Evans 2008). Given some reasonable assumptions about 
the temperature and morphology of a thermal shell, comprehensive 
radio spectral coverage as the ejecta transition from being optically 
thick to optically thin enables observers to infer the density pro¬ 
file of the nova shell and estimate the total mass ejected. Resolved 
images in radio allow for further constraints on the morphology, 
temperature, and evolution of the thermal shell. However, only a 
handful of classical novae have been well observed at multiple ra¬ 
dio frequencies from the optically thick rise to the optically thin 
decay, including HR Del (Hjellming 1979), V1500 Cyg (Seaquist 
et al. 1980), and V1974 Cyg (Hjellming 1996). This lack is partic¬ 
ularly noticeable in the first few months immediately following the 
outburst, when the thermal radio emission is expected to be faint. 

Technological improvements to the Karl G. Jansky Very Large 
Array (VLA) and other radio telescopes within the past five years 
have vastly increased the sensitivity of radio observations. Our 
group has been using the updated capabilities of the VLA in con¬ 
cert with optical and X-ray observations to create a more complete 
and multi-wavelength picture of the mass-loss history of novae. 
V1723 Aql was the first classical nova to be discovered in outburst 
after the upgrade to that facility. Beginning two weeks after the ini¬ 
tial discovery of the eruption of V1723 Aql in 2010 September, the 
observations of Krauss et al. (2011), in combination with the addi¬ 
tional radio monitoring presented here, have produced one of the 
best, most detailed radio light curves of a classical nova to date. 

Discovered in outburst at its optical maximum (mag 12.4) on 
2010 September 11, and with a magnitude of fainter than 14.0 on 
2010 September 10, we take that discovery of V1723 Aql as the 
start of the outburst, to = MJD 55450.5 (Yamanaka et al. 2010; 
Balam et al. 2010). Optical observations revealed V1723 Aql to be 
a typical ‘fast nova’ (Gaposchkin 1957) of the Fe II spectral type, 
fading by two magnitudes within 20 days (Yamanaka et al. 2010). 
Despite the lack of extensive optical or IR followup, the observa¬ 
tions in the first few weeks of the erruption showed rapid reddening, 
implying the formation of dust particles as early as 20 days after the 
start of the outburst (Nagashima et al. 2013). Emission lines of Hir 
on the day of discovery had full width at zero intensity of 3000 
km s“*, implying a maximum expansion velocity Vobs =1500 km 
s“' (Yamanaka et al. 2010). The nova was neither embedded in the 
wind of an M-giant mass donor nor a recurrent nova - making it 
representative of the majority of novae (Nagashima et al. 2013). 

Although V1723 Aql showed no unusual aspects in its early 
optical behaviour, the VLA observations revealed an unexpected 
radio flare peaking roughly 45 days after the start of the optical 
outburst (Krauss et al. 2011). Krauss et al. (2011) found that dur¬ 
ing this flare, the radio flux density rose as t^'^, far more rapidly 
than the t^ expected from freely expanding, isothermal, optically 
thick thermal ejecta. The radio spectrum during the flare was also 
inconsistent with that of an optically thick thermal source of a fixed 
temperature (Krauss et al. 2011). With only 175 days of radio moni¬ 
toring, Krauss et al. (2011) were not able to conclusively determine 
the origin of this flare, or whether V1723 Aql would evolve accord¬ 
ing to the more standard picture of radio emission from a classical 
nova after this unexpected activity. We label the observations be¬ 
tween the start of the eruption and day 175, presented in Krauss et 
al. (2011), as the early-time observations, and the new observations 


presented here (between day 198 and day 1281) as the late-time 
observations. 

We describe below how the combination of the early-time and 
late-time radio observations reveals that whereas the late-time ob¬ 
servations are dominated by the expected thermal emission from a 
freely expanding nova shell, there is strong evidence that the early- 
time flare was dominated by non-thermal emission. In section 2 we 
present our late-time observations and data reduction procedure; 
in section 3, we describe the observational findings and data anal¬ 
ysis; in section 4, we introduce a thermal model for the remnant 
after day 200, explore the possibility of the early-time emission be¬ 
ing produced by relativistic particles accelerated in internal shocks, 
and propose that the ejecta from V1723 Aql had a morpho-kinetic 
structure with internal shocks like that of y-ray producing nova 
V959 Mon. In section 5, we highlight our conclusions. 


2 OBSERVATIONS AND DATA PROCESSING 

Building on the work of Krauss et al. (2011), who used 
VLA observations during the first six months of the outburst of 
V1723 Aql to identify the bright early-time flare, we present our 
observations of V1723 Aql with the VLA between 2011 March 
(day 198) and 2014 March (day 1281) - from six months to three 
and a half years after the start of the eruption (see Table 1). For 
completeness. Table 1 additionally lists the early-time observations 
of Krauss et al. (2011). The late-time observations were taken using 
the VLA WIDAR correlator at a range of frequencies spanning 2.5- 
36.5 GHz with bandwidths of 2048 MHz, distributed in two tune¬ 
able 1048-MHz sidebands and using 8-bit samplers. The central 
frequencies in each band may have had some slight variation due 
to frequency dependent flagging in post-processing to minimize ra¬ 
dio frequency interference (RFI). At each frequency, on-source ob¬ 
servations were alternated with observations of a nearby phase ref¬ 
erence calibrator (J1851-^0035, J1832-1035, or J1822-0938), with 
an additional observation of a standard flux calibration source (ei¬ 
ther 3c286 or 3c48). We used observations of the unpolarized cal¬ 
ibrator J0319-I-4130 (3c84) for polarization leakage solutions. The 
late-time data presented here comprise 29.5 hours of observations 
(14.9 hours on source V1723 Aql) and 658 GB of raw data over 
three years and two full cycles of VLA configurations. Including 
the earlier programs presented in Krauss et al. (2011), the data to¬ 
tal 1604 GB and 59.5 hours (28.2 hours on source) between 2010 
September 25 and 2014 March 15. 

We processed all data using NRAO’s Common Astronomi¬ 
cal Processing System (CASA) or Astronomical Image Processing 
System (AIPS). As is standard procedure, we obtained a bandpass 
solution and used it to solve for gain amplitude and phase, then ap¬ 
plied these solutions to the source, using the ‘Perley-Butler 2010’ 
flux density scale. Imaging was completed using the CLEAN algo¬ 
rithm (Hogbom 1974) with Briggs weighting, with the additional 
use of superresolution techniques (Staveley-Smith et al. 1993) on 
day 1281 to resolve small structures during that observation. Radio 
flux densities were measured by fitting two dimensional gaussians 
to the emitting region in the image plane with the task imfit. We 
estimated uncertainties in flux density by adding in quadrature the 
error from the gaussian fit and systematic errors of 1% at frequen¬ 
cies below 19 GHz and 3% at those above, with the results pre¬ 
sented in Table 2 and Figs 1 and 2. All error bars correspond to Icr 
uncertainties. 
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Table 1. VLA observations of V1723 Aql, between 2010 September and 2014 March. 


Observation Date 

MJD 

Day" 

Epoch 

VLA Program ID 

Configuration 

Observed Bands 

Time* (min) 

Early time observations (previously presented in 

Krauss et al. 2011) 




2010 September 25 

55464.1 

13.6 

0 

lOB-200 

DnC 

C, Ka 

57 

2010 October 3 

55472.1 

21.6 

1 

lOB-200 

DnC 

C, Ka 

21 

2010 October 6 

55475.1 

24.6 

2 

lOB-200 

DnC/C 

C, Ka 

58 

2010 October 15 

55483.1 

32.6 

3 

lOB-200 

C 

C, Ka 

48 

2010 October 18 

55487.1 

36.6 

4 

lOB-200 

c 

L, C,X 

43 

2010 October 19 

55488.1 

37.6 

4 

lOB-200 

c 

C, Ka 

28 

2010 October 24 

55493.1 

42.6 

5 

lOB-200 

c 

L, C,X 

43 

2010 October 29 

55497.0 

46.5 

6 

lOB-200 

c 

L, C, X, Ka 

71 

2010 November 7 

55507.8 

57.2 

7 

lOB-200 

c 

L, C, X, Ka 

74 

2010 December 3 

55533.9 

83.4 

8 

lOB-200 

c 

L, C, X, Ka 

84 

2010 December 18 

55548.9 

98.4 

9 

lOB-200 

c 

L, C,X 

52 

2010 December 21 

55551.7 

101.2 

9 

lOB-200 

c 

C, Ka 

34 

2011 January 14 

55575.8 

125.3 

10 

lOB-200 

c 

L,C,X, Ka 

70 

2011 January 30 

55591.8 

141.3 

11 

lOB-200 

CnB 

L,C,X, Ka 

63 

2011 March 5 

55625.5 

175.0 

12 

llA-254 

B 

C, X, Ka 

51 

Late time observations 

2011 March 28 

55648.6 

198.1 

13 

llA-254 

B 

C, X, Ka 

52 

2011 August 2 

55775.2 

324.7 

14 

llA-254 

A 

C, X, Ku, K, Ka 

116 

2011 August 27 

55800.2 

349.7 

15 

llA-254 

A 

C, X, Ku, K, Ka 

72 

2011 November 11 

55876.1 

425.6 

16 

llB-170 

D 

C, X, Ku, K, Ka 

58 

2012 January 14 

55940.8 

490.3 

17 

llB-170 

DnC 

C, X, Ku, K, Ka 

58 

2012 April 19 

56037.0 

586.5 

18 

llB-170 

C 

C, X, Ku, K, Ka 

58 

2012 December 8 

56269.0 

818.5 

19 

12B-226 

A 

C, X, Ku, K 

20 

2012 December 22 

56283.6 

832.5 

20 

12B-226 

A 

X, Ka 

45 

2013 March 27 

56378.4 

927.9 

21 

13A-465 

D 

C, X, Ku, K, Ka 

106 

2013 August 30 

56534.1 

1083.6 

22 

13A-461 

C 

L, C, X, Ku, K, Ka 

98 

2013 November 30 

56626.7 

1176.5 

23 

13A-465 

B 

L, S, X, Ka 

101 

2014 March 15 

56731.5 

1281.0 

24 

13A-465 

A 

X, Ku 

112 


“ Days after initial detection of outburst, 2010 September 11. 
* Total time of observations on-source. 


3 DATA ANALYSIS AND RESULTS 

3.1 Flux density and spectral evolution 

Around day 45, the flux densities reached a local maximum 
and subsequently declined at all frequencies, with another local 
maximum at all frequencies between days 200 and 400 (see Table 2 
and Fig. 1). We refer to the radio brightening that peaked around 
day 45 as the early-time flare. There was an additional local max¬ 
imum between these two bright peaks at the higher frequencies. 
The flux densities at the highest frequencies reached higher values 
than those at lower frequencies during both major peaks. For ex¬ 
ample, during the final local maxima (between days 200 and 400), 
the 36.5 GHz flux density peaked at 12.80 ± 0.39 mjy, whereas the 

5.2 GHz flux density peaked at 1.14 ± 0.05 mJy. Moreover, the flux 
densities at Aa-band (28.5/36.5 GHz) reached their highest values 
on day 198, whereas the flux densities at lower frequencies peaked 
later, with C-band flux densities (5.2/6.8 GHz) peaking on day 426 
or 490. 

During the late-time period, the radio spectrum evolved from 
one that rose with frequency to one that fell slightly with frequency 
(see Fig. 2 and Table 3). To obtain spectral indices and la un¬ 
certainty estimates for those spectral indices for each observation 
(5v cx v“, where 5v is the flux density at frequency v, and a is the 
spectral index), we fit data with a power-law using IDL’s Tinfit’ 
function, which minimizes the statistic, with the data and er¬ 
ror bars in log space. With the simple single power-law model, we 


found acceptable fits ixl <1-5, where is the reduced value) 
only for days 198 and 818, between which time the spectral index 
fell from a = 1.47 ± 0.05 to a = -0.22 ± 0.08. 

The late-time transition from a rising to a falling spectrum oc¬ 
curred via the appearances of first one, and then a second, spec¬ 
tral break. We refer to the three distinct portions of the spectrum 
when two breaks were present as the low-frequency, transition, and 
high-frequency parts of the spectrum. We refer to the corresponding 
spectral indices as aiow, tttrans, and a^, and the frequencies at which 
the breaks between the regions of different spectral index occur as 
vit (for the break between the low-frequency and transition part of 
the spectrum) and Vth (for the break between the transition and high- 
frequency portion of the spectrum). For epochs where both a poor 
X^ fit and the residuals of the simple power-law model suggested a 
spectral break, we fit the data with broken power law models with 
one or two breaks. Fig. 2 shows the spectra and fits, and Table 3 
reports the power law models that combined to minimize with 
the break frequencies indicating where the power-law models in¬ 
tersect. The uncertainties for the break frequencies correspond to 
the range of frequencies at which the power laws intersected given 
a ±lcr range for the spectral indices. We were able to find accept¬ 
able fits for a single break power-law model on day 325, and for a 
double-break power law model on days 426 and 587. 

Fig. 2 and Table 3 show the evolution of the spectrum and its 
component segments. The low-frequency spectral index aio„ typ¬ 
ically had values of approximately 1.4, tttrans typically had values 
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Table 2. Flux densities for late-time observations of V1723 Aql 


Epoch 

Day 

2.5 GHz (S band) 

3.3 GHz (S band) 

Observed flux density (mjy) 

5.2 GHz (C band) 6.8 GHz (C band) 

8.7 GHz (X band) 

9.5 GHz (X band) 

13 

198 

* 


0.75 ± 0.02 


1.53 ±0.03 


14 

325 



0.89 ± 0.03 

1.32 ±0.03 

1.88 ±0.03 


15 

250 



0.91 ± 0.03 

1.33 ± 0.03 

1.97 ±0.03 


16 

425 



1.14 ± 0.05 

1.59 ±0.04 

2.09 ± 0.04 


17 

490 



1.10 ± 0.03 

1.60 ±0.05 

1.79 ±0.03 


18 

587 



0.92 ± 0.02 

1.30 ±0.02 

1.39 ±0.03 


19 

818 



0.70 ± 0.05 

0.76 ± 0.06 



21 

928 



1.09 ±0.10 

0.61 ± 0.03 

0.77 ± 0.04 

0.71 ± 0.02 

22 

1084 



0.43 ± 0.02 

0.48 ± 0.02 

0.44 ± 0.02 

0.45 ± 0.02 

23 

1177 

0.42 ± 0.04 

0.25 ± 0.02 



0.41 ± 0.05 

0.34 ± 0.01 

Epoch 

Day 



Observed flux density (mJy) 





12.9 GHz (Ku band) 

16.1 GHz (Ku band) 

20.1 GHz (K band) 

25.6 GHz (K band) 

28.5 GHz (Ka band) 

36.5 GHz (Ka band) 

13 

198 



5.48 ±0.17 

7.76 ± 0.23 

9.02 ± 0.28 

12.81 ±0.39 

14 

325 

3.20 ±0.12 

3.94 ±0.10 

4.87 ±0.16 

5.35 ± 0.17 

5.84 ±0.18 

6.76 ± 0.22 

15 

350 

2.98 ± 0.05 

3.59 ± 0.05 

4.20 ±0.13 

4.30 ±0.14 

4.43 ±0.14 

5.38 ± 0.20 

16 

426 

2.42 ± 0.03 

2.71 ±0.08 

3.03 ±0.10 

3.03 ± 0.10 

2.87 ± 0.09 

3.12 ±0.11 

17 

490 

2.10 ±0.03 

2.07 ± 0.03 



2.42 ± 0.09 

2.05 ± 0.10 

18 

587 

1.53 ±0.04 


1.87 ±0.08 

1.72 ±0.07 

1.85 ±0.08 

1.69 ±0.10 

19 

818 

0.66 ± 0.06 

0.56 ± 0.05 


0.43 ± 0.09 



20 

832 





0.71 ±0.13 

0.62 ±0.19 

21 

928 

0.53 ± 0.02 

0.59 ± 0.02 

0.54 ± 0.03 

0.56 ± 0.03 

0.55 ± 0.02 

0.54 ± 0.03“ 

22 

1084 

0.43 ± 0.03 

0.42 ± 0.02 

0.51 ±0.08 

0.32 ± 0.03 

0.38 ± 0.05 

0.35 ± 0.04“ 

23 

1177 





0.23 ± 0.02 

0.19 ±0.02“ 

24 

1281 

0.27 ± 0.01“ 

0.26 ± 0.02* 






* indicates no observation at this frequency on this day. 
" Values are reported for 13.5 GHz. 

^ Values are reported for 14.5 GHz. 

^ Values are reported for 32.5 GHz. 


Table 3. Spectral indices. For days where we fit models with spectral breaks, vit refers to the break between the low-frequency and 
transition parts of the spectrum, and Vl\^ refers to the break between the transition and high-frequency portions of the spectrum. 
aiow» ti'trans. and cTh are the corresponding spectral indices for these segments. shows the density profile n oc l/r^ associated 
with the spectral index a for a spherical shell transitioning between optically thick and optically thin thermal emission (Wright 
& Barlow 1975), and;^^ and df are the reduced chi-square statistic and the number of degrees of freedom, respectively. 


Epoch 

Day 

model* 

U'/ow 

^trans 


VI, ( GHz ) 

v,h ( GHz ) 


df 

P 

13 

198 

SPL 

1.48 ± 0.02 





0.73 

4 


14 

325 

BPL 

1.44 ± 0.05 

0.64 ± 0.05 


13.6 


0.45 

1 

2.04 ± 0.06 

15 

350 

BPL 

1.54 ± 0.06 

0.52 ± 0.03 


11.3 


2.75 

7 

1.92 ± 0.03 

16 

426 

DBPL 

1.23 ± 0.19 

0.41 ± 0.04 

0.03 ± 0.08 

8.3 

21.2 

0.60 

7 

1.84 ± 0.03 

17 

490 

DBPL 

1.37 ± 0.15 

0.24 ± 0.04 

- 0.64 ± 0.23 

7.3 

27.8 

2.53 

5 

1.69 ± 0.03 

18 

587 

DBPL 

1.26 ± 0.10 

0.30 ± 0.05 

- 0.08 ± 0.19 

6.7 

21.2 

1.00 

6 

1.74 ± 0.04 

19 

818 

SPL 



- 0.22 ± 0.08 



1.26 

3 


20 

832 

SPL 



- 0.47 ± 1.34 



- 

0 


21 

928 

SPL 



- 0.19 ± 0.02 



8.40 

8 


22 

1084 

SPL 



- 0.10 ± 0.04 



2.12 

8 


23 

1177 

SPL 



- 0.08 ± 0.03 



24.19 

4 


24 

1281 

SPL 



- 0.22 ± 0.77 



- 

0 



*SPL - single power law fit 

BPL - broken power law fit (one break) 

DBPL - double broken power law fit (two breaks) 


between 0.2 and 0.6, and typically had values of approximately 
-0.2. The spectrum on day 490 (see Fig. 2) shows the character¬ 
istic low-frequency, transition, and high-frequency spectral regions 
most clearly. On day 198, the spectrum was well described by a 
simple power-law with a low-frequency type spectral index. On day 
325, the lowest frequencies maintained a steep rising spectrum of 


triow = 1.44+0.05; above 16.1 GHz the spectrum became shallower, 
with t^trans = 0.64 ± 0.05. The spectra continued to show this break 
between (3fiow and o^trans in observations between days 325 and 587, 
during which time v\i decreased from 13.6 GHz to 6.7 GHz. Ob¬ 
servations between days 426 and 587 showed the additional break 
at higher frequencies, Vjh, above which the spectrum flattened to 
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Figure 1. Observed flux densities of V1723 Aql between 2010 September and 2014 March. Error bars are as reported in Table 2, but may be too small to be 
visible. 



ffh = 0.13 ± 0.13 on day 426, orh = -0.64 ± 0.23 on day 490, and 
tth = -0.08 ± 0.19 on day 587. While the very small uncertainties 
in our flux densities reveal some deviation from a single power law 
during the last five observations, the spectrum was nevertheless ap¬ 
proximately described by a fairly flat single power law during that 
time. 

3.2 Resolved images 

V1723 Aql was resolved at the highest observed radio fre¬ 
quency on day 832 (Ka-band) and on day 1281 (Ku-band) (see 
Fig. 3, panel a). The remnant was therefore large enough to be spa¬ 
tially resolved, with the difference between the integrated and the 
peak flux density greater than five times the root mean square error 
(RMS). Fig. 3 shows images from the second and third of the three 
observations during which the VLA was in A configuration (which 
gives the highest spatial resolution). The image from day 832 (2012 
December 22) at 32.4 GFlz shows that the radio emitting material 
was elongated in the NW/SE direction and had visible substruc¬ 
ture. The full width half maximum of the brightness profile along 
the major axis of the fitted gaussian was 160 ± 10 milliarcseconds 
(mas), with the FWFIM along the minor axis of 110 ± 20 mas. 
The position angle (PA) of the major axis wasl60° ± 10° (where PA 
is the east of north). The clean beam size was 90 mas by 60 mas. 


and the peak flux density was 130 ± 20 yujy/beam, with an RMS of 
22.6 /jJy/beam. To more accurately estimate the size scale of the 
source, we fit the brightness profile, deconvolved from the synthe¬ 
sized clean beam, to a Gaussian. Taking the width of the Gaussian 
along the major axis at a flux density of twice the RMS as an in¬ 
dicator of size, the size scale along the major axis on day 832 was 
250 ±10 mas. The substructure showed two resolved components 
with brightness peaks separated by 70 mas. They were aligned per¬ 
pendicular to the major axis of the elongated material, with a PA of 
about 80°. 

Taking into account the expected expansion by 50% between 
day 832 and day 1281, the size of the radio-emitting remnant on day 
1281 indicates that the material that was most extended on day 832 
had faded, and only the denser inner structure remained detectable. 
By day 1281 (2014 March 14; using robust = 0.5 in the CLEAN 
algorithm), the resolved material was elongated in the NE/SW di¬ 
rection, instead of the NW/SE elongation seen on day 832. The 
fitted gaussian had a FWFIM of 180 ± 10 mas along the major axis 
(which had a PA of 35° ± 5°), and a FWFIM along the minor axis of 
133 ± 7 mas. The peak flux density was 118.5 ± 3.5 /rjy/beam, with 
RMS of 3.32 yujy/beam, resulting in a size scale of 408 ± 2 mas. 
With a clean beam size of 160 mas by 130 mas, no substructure 
was resolved. If we consider the two components visible on day 
832, and extrapolate from their 70 mas separation, we would have 
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Figure 2. Spectra of V1723 Aql between 2011 March and 2014 March, with lines showing the best-fit spectral indices a (where Sy oc v“). For epochs where the 
spectral index appears to be a function of frequency, spectral indices are listed from lowest frequency range to highest. See Table 3 
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expected a separation of 108 mas on day 1281 - which was not 
within the formal resolution limits of this observation. Therefore, 
we additionally used superresolution techniques (Staveley-Smith 
et al. 1993) to examine the more detailed structure of the source 
on this day, deconvolving the source with the slightly smaller beam 
size of 90 by 90 mas. Using this technique, we uncovered a ring¬ 
like structure with dimensions ~ 0.22 x 1.5 arcsec and PA 45° (see 
Fig. 3, panel b), with an RMS of 7.0 yujy/beam. 


4 DISCUSSION 

We propose that the radio emission from V1723 Aql - when 
considered in the context of published optical and X-ray observa¬ 
tions - is most consistent with a picture in which the nova initially 
generated a slow outflow and then a faster flow. We argue below 
that the slow flow had a speed of a few hundred km s“*, and that it 
was likely concentrated toward the equatorial plane of the binary. 
The faster flow had a maximum speed of at least 1500 km s“*, 
and it collided with the early, slow ejecta within a few weeks of 
the start of the eruption, producing the early-time flare. The colli¬ 
sions led to shocks that heated some of the ejecta to temperatures 
above 10^ K, and also accelerated particles to relativistic speeds. 
The slow-moving material probably collimated the faster flow into 
a biconical structure. Fig. 4 shows a cartoon of this toy model. The 
scenario that we propose for V1723 Aql is consistent with the one 


that Chomiuk et al. (2014) suggested for the y-ray producing nova 
V959 Mon. 

In the discussion below, we present evidence that supports the 
picture of a fast flow that collided with, and was shaped by, a pre¬ 
ceding slow outflow. To diagnose the flows without the complexity 
of non-thermal emission from the early-time shock, we begin in 
section 4.1 by describing how the brightness temperature and radio 
spectral evolution at late times were consistent with a freely ex¬ 
panding thermal plasma with a temperature of around 10'* K. The 
density profile and radio images of this radio-emitting plasma sug¬ 
gest that it had a roughly bipolar or biconical outer structure, and a 
shell-like inner structure with a likely equatorial density enhance¬ 
ment. In section 4.2, we show that the early-time radio flare was 
probably dominated by synchrotron emission. Given the low X-ray 
flux during the flare (Krauss et al. 2011) and high brightness tem¬ 
peratures, it cannot be attributed to thermal emission alone. In sec¬ 
tion 4.3, we discuss how V1723 Aql relates to other classical novae 
in its behaviour and shape, including those that produce y-rays. 

4.1 Origin of the Late-Time Radio Emission: Expanding 
thermal shell 

4.1.1 General physical constraints 

We obtained a first estimate of the distance to V1723 Aql by 
using our maximum extension of the remnant on day 832, when 
our image had the highest spatial resolution (smallest clean beam). 
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Table 4. Total ejecta mass M and minimum velocity (v/) as determined using the spectral breaks visible on 
days 426, 490, and 587. Mass estimates vary based on the density profile assumed. 


Epoch 

Day 

V/ (km s ^) 

C = V,7Vobs 

M (Mq), ncc 1/r^ 

M {Mq), n cc l/r“ 

M (Mq), ncc \jp 

16 

426 

360 ± 20 

0.24 ± 0.01 

2.7±0.3xl0-"‘ 

1.6±0.2xl0^'‘ 

0.9 ±0.2x10-'* 

17 

490 

210 ± 10 

0.14 ±0.01 

4.0±0.4xl0-"‘ 

1.6±0.2xl0^'‘ 

0.6 ± 0.1 xlO-”* 

18 

587 

190 ± 20 

0.13 ±0.01 

4.6 ± 1.6 xlO-'^ 

1.7 ±0.7x10^“* 

0.7 ±0.3x10-'* 


Figure 3. Top panel: Spatially resolved image of V1723 Aql taken on 2012 De¬ 
cember 22 in Ka-hiccvd (32.4 GHz). When fitted to a two-dimensional gaussian, 
the brightness profile has a FWHM along the major axis of 160 ± 10 milliarcsec- 
onds (mas) and along the minor axis of 110 ± 20 mas, where the major axis has a 
position angle of 160° ± 10° east of north. Contour levels are at [-2,2,3,4,5]xRMS 
with RMS= 2.26 x 10“^ Jy. There are two resolved components roughly perpen¬ 
dicular to the major axis of the gaussian profile, with a separation of 70 mas at 
an angle of 82° E of N. Bottom panel: Spatially resolved image of V1723 Aql 
taken on 2014 March 14 in ^w-band (14 GHz). When deconvolved with a beam 
size of 90 X 90 mas, we see a ring-like structure. Contour levels are at [-2, 2, 
2 V2, 4, 4 V2, 8, 8 V2, 16, 16 V2, 32]xRMS with RMS=7.01 x 10“^ Jy. 
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Since by day 832, some of the outermost material had become op¬ 
tically thin and therefore the remnant might not have a sharp outer 
edge, we took the maximum size scale to be the extent along the 
major axis of the gaussian fit, down to a flux density of twice 
the RMS. When we compare the size scale of 251 ±7 mas with 
the expectation for ejecta expanding with a maximum velocity of 
Vobs = 1500 km s“' (Yamanaka et al. 2010), we obtain a distance of 
5.7 ± 0.4 kpc, or ~6 kpc. The assumption that the maximum expan¬ 
sion speed in the plane of the sky is equal to the maximum radial 
velocity is valid if the remnant is approximately spherical or, for a 


Figure 4. We propose a picture where a fast, bipolar outflow is shaped by a 
slow moving equatorial torus. Shocks arose in the collision region between 
the slow and fast flows. 



bipolar flow, has wide opening angles (as is suggested by the 1 jd 
density profile, see below). 

Other methods of distance determination using optical data are 
either not possible or extremely unreliable. There are two means 
of nova distance determination using optical data. The most com¬ 
mon utilized method is the Maximum Magnitude Rate of Decline 
(MMRD) method (Della Valle & Livio 1995); the other method 
involves utilizing high-resolution spectra to identify absorption 
lines from intervening interstellar clouds of known distance (e.g. 
Chomiuk et al. 2012). No high resolution spectra of V1723 Aql 
were available, making it impossible for us to use the latter method. 
The MMRD technique has been the subject of serious critiques in 
recent publications (e.g. Kasliwal et al. 2011; Cao et al. 2012). 
Further, the episode of dust formation observed in V1723 Aql 
confounds the determination of t 2 , rendering the already dubious 
MMRD method useless for distance determination. 

Using d~6 kpc, we find that the brightness temperature after 
day 200 was consistent with the late-time radio emission emanating 
from a thermal remnant with T ~ 10"* K, where T is the electron 
temperature. Brightness temperature is 




Sy(t)dD^ 
2;rfci,v2(Vobsf)^ ’ 


( 1 ) 


where c is the speed of light, is Boltzmann’s constant, t is the 
time since eruption, and v^bs is the observed maximum radial ve¬ 
locity (e.g. Bode & Evans 2008). Using our values from day 200, 
for example, we find that T\, = 1.5 ± 0.4 x 10“* K at 5.2 GHz, 
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which is consistent with ~ 10"^ K. Ionized ejecta from novae are typ¬ 
ically found to have electron temperatures of around 10"* K, which 
is also approximately the equilibrium temperature for photoionized 
plasma that is cooled by forbidden line emission (Seaquist & Pali- 
maka 1977; Taylor et al. 1987; Bode & Evans 2008). This ~10‘* K 
ejecta temperature can be sustained by photoionization heating by 
the WD for up to a year after the nova eruption (Cunningham, 
Wolf & Bildsten 2015). A brightness temperature around this value 
therefore provides support for thermal bremsstrahlung as the dom¬ 
inant emission mechanism at late times in V1723 Aql. 

Moreover, the detailed evolution of the radio spectrum be¬ 
tween days 198 and 350 (approximately 6 to 12 months after the 
start of the outburst) places constraints on the density profile in the 
outermost regions of the expanding ejecta. We consider the point 
at which the free-free optical depth at a given frequency (t,) is 
unity, where the emitting material transitions from optically thick 
to optically thin. The high spectral index of aiow ~ 1-5 on day 198 
indicates that at that time, the ejecta were dense enough to be op¬ 
tically thick for frequencies of up to at least 36.5 GHz. Therefore, 
these Ty = \ surfaces were very close to the outer boundary of the 
ejecta for all frequencies at which we observed on day 198. The 
appearance of the first spectral break on day 325, however, reveals 
that by that time, the density of the ejecta had dropped enough that 
at least part of the ejecta were beginning to become optically thin 
for frequencies above vu. Thus, the radii of the Ty = \ surfaces at 
frequencies above Vj, were no longer increasing as fast as that of the 
ejecta. The spectral index above the first spectral break (fftrans) is di¬ 
rectly related to the density profile in these transitioning ejecta. For 
an isothermal shell with a density profile of the form n{r) oc 1 /r^ 
(where n is density and r is distance from the central binary), the 
spectral index is given by tttrans ~ for yS > 1.5 (Wright & 

Barlow 1975). For a measured tttrans, we thus have yS = 2o‘-iOtttr»s 
(see Table 3). Our measurement of ttirans ~ 0.6 on days 325 and 
350, when only the outermost portions of the ejecta were begin¬ 
ning to become optically thin, suggests that the density profile in 
the outermost part of the ejecta was proportional to 1/r^. 

The appearance of a second break in the spectrum on day 426, 
roughly 14 months after the start of the eruption, shows that the 
ejecta had become optically thin at frequencies above v,),. Optically 
thin thermal emission has a spectral index of a = -0.1 (e.g. Bode 
& Evans 2008),as we approximately found above Vth- Between the 
two breaks at V;, and v,;,, the spectral index of the transitioning re¬ 
gion tttrans provides information about the density profile deeper 
within the ejecta. With a value of approximately 0.3 on days 490 
and 587, when these two breaks are clear, rttrans suggests an average 
density profile deep within the remnant of roughly l/r''^ - some¬ 
what flatter than in the outermost regions of the remnant. 

The observed evolution of rttrans is unlikely to be due to a tem¬ 
perature gradient. Although a temperature gradient can in princi¬ 
ple affect the spectral index, the effect drops out for a 1 /r^ density 
profile (Cassinelli & Hartmann 1977; Barlow 1979). Moreover, a 
radial decrease in the ionized gas fraction would also only serve to 
increase the steepness of the spectrum during the transition period 
(Wright & Barlow 1975), in contrast to what we observed. There¬ 
fore, the observed flattening of the rttram in the interior of the ejecta 
is unlikely to have been due to a temperature or ionization gradient. 

For a more promising avenue to explain the declining fftransi 
Wright & Barlow (1975) note that non-spherical geometries can 
cause the total number density of the gas to fall off less rapidly than 
n oc for a constant, finite duration flow. For example, a disk 
geometry would fall off as « cx and in a cylindrical geometry 


n would be independent of r (Wright & Barlow 1975). Therefore, 
that for V1723, Aql trirans appears to flatten slightly with time may 
indicate that the inner remnant is more confined towards the equa¬ 
torial plane than the outer remnant - as in the case of the biconical 
flow collimated by a dense, slow torus in V959 Mon (Chomiuk et 
al. 2014). If this is the case, a bi-polar outflow could result in an 
over-estimate of mass by up to a factor of two (Ribeiro et al. 2014). 
Nevertheless, since the density profile of 1/r^ in the outer regions of 
the ejecta suggest that any bipolar flow must have a wide opening 
angle, we can nevertheless use a spherical model to obtain reason¬ 
able estimates of the ejecta parameters. 

The flux density near v* furnishes a rough estimate of the in¬ 
ner radius of the ejecta on the dates when this break frequency was 
observable. The ejecta are optically thin, with optical depth Ty < 1, 
at frequencies above v*. They are optically thick, with Ty > 1 at fre¬ 
quencies below Vth- Integrating the absorption coefficient for radio 
emission at v* along the line of sight through the ejected envelope 
therefore gives exactly Ty = 1, suggesting that the photosphere lies 
near the inner edge of the ejected shell for Vth. At the frequency Vthi 
the rough angular size of the photosphere is 


e ^ / Sy n 

arcsec ^mjyj ^1200Kj ^GHz 


( 2 ) 


(Rybicki & Lightman 1979), where 9 is the full angular diame¬ 
ter of the photosphere. Using the break frequencies Vth and flux 
densities at these frequencies from days 426, 490, and 587 (see 
Fig. 2), and taking T* » 10"* K, we find that the inner edge of 
the ejecta had an approximate radius of between 10 and 15 mas 
during these observations, which corresponds to an expansion 
speed V, of the slowest, innermost ejecta of between 200 and 
400 km s“*(d/6 kpcjiTj/lO"* K),“'^^, assuming all the ejecta were 
homologous and ejected at tg. Although these estimates do not pro¬ 
vide any evidence that the inner radius increased between days 426 
and 587, the estimates are uncertain enough that they do not place 
meaningful constraints on the actual evolution of the inner radius. 
The estimates do, however, indicate that the radio-emitting ejecta 
were geometrically thick (given that the outermost material moved 
away from the central binary at approximately 1500 km s“'), 
though not uniquely or unusually so (e.g. Seaquist & Palimaka 
1977). 

The radio spectra during the second year of the eruption also 
enable us to make an estimate of the ejecta mass. We can obtain the 
emission measure of the ejecta at the Ty = 1 radio photosphere at 
Vth from the expression 


Ty = 8.235 X 10^^ 



/ \ 
\ cm-® pc / 


(3) 


where Ty is the free-free optical depth at frequency v and EM^ad = 
” rp-^r) dr is the emission measure integrated from the inner ra¬ 
dius Ri to the outer radius Rg (Lang 1980). With the emission mea¬ 
sure, inner and outer radii, and density profile in hand for the ob¬ 
servation with the clearest break frequency v,h, we can estimate the 
mass of the ejecta. For the purpose of this mass estimate, we take a 
spherically symmetric shell and density proportional to 1/r^, as es¬ 
timated from iTtrans in the outer regions of the ejecta. Using the inner 
radius on day 426 from above, which gives the velocity of the inner 
edge as v, = 360 km s“' and outer radius assuming a maximum 
expansion velocity of 1500 km s^', we find an approximate ejecta 
mass of 1.6(±0.2) x 10“"^ Mq, which is of the same order of mag¬ 
nitude as is typically found in nova ejecta mass estimates (Bode & 
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Evans 2008). Repeating this process for other days on which mul¬ 
tiple breaks were visible results in similar mass estimates - on day 
490 we find a mass of 1.6(±0.2) x 10“** Mq, and on day 587 we find 
a mass of 1.7(±0.7) x 10“"* Mq (see Table 4). 

4.1.2 The Hubble Flow model 

The gross properties of the radio-emitting ejecta at late times 
(after day 200) were consistent with expectations from the so- 
called “Hubble-Flow” model (Seaquist c& Palimaka 1977; Hjellm- 
ing 1979). In this model, a spherically symmetric shell is ejected 
from the WD instantaneously, but with a range of velocities, so 
that the velocity is proportional to the distance from the WD. A 
1/r^ density profile, as we inferred from atrans, corresponds to the 
mass being distributed uniformly over the velocities. Such a model 
has been applied to many historical novae, including VI974 Cyg 
(Hjellming 1996), V1500 Cyg (Seaquist et al. 1980), QU Vul (Tay¬ 
lor et al. 1988), and V723 Cas (Heywood et al. 2005). Given the 
brightness temperature after day 200 of near 10"* K, and the evolu¬ 
tion of the radio spectrum from initially rising with frequency (as 
expected for an optically thick shell) to falling slightly with fre¬ 
quency (as expected for an optically thin shell) via the appearance 
of two distinct spectral breaks, we conclude that it is reasonable to 
use a Hubble-flow model to parametrize the ejecta that produced 
the late-time radio emission from V1723 Aql. 


Increasing the total ejecta mass in the model also causes an increase 
in flux density at all times, and causes a slower, later transition from 
optically thick to optically thin emission at all frequencies. 

While our observations are generally consistent with anoc 
density profile, we can additionally demonstrate the sensitivity of 
the total ejecta mass to the distribution of mass within the nova 
shell. For example, if we again use the observed spectral breaks 
to estimate the inner and outer radii of the shell and calculate the 
emission measure, we can estimate the mass of the shell with dif¬ 
ferent density profiles. For example, by using the spectral breaks 
on day 587, we estimate an ejecta mass of M ~ 4.6 x 10 “"^Mq 
and M ~ 0.7 x 10“‘'Mo for densities profiles of 1/r and 1/d, re¬ 
spectively (see Table 4). If the ejecta from VI723 Aql were clumpy 
rather than smooth, the total ejecta mass would be reduced. Clump¬ 
ing increases the flux density of an isothermal shell by up to a factor 
of where / is the volume filling factor (Abbott et al. 1981). 
However, clumping would not influence the spectral index, nor 
the overall shape of the radio light curve, as uniformly distributed 
clumping would effect the free-free flux similarly at all frequencies 
(Abbott et al. 1981; Leitherer & Robert 1991; Nelson et al. 2014). 


4.2 Origin of the early time flare: particles accelerated in 
shocks within the ejecta 

4.2.1 Brightness temperature of the ejecta 


Although the Hubble-flow model does not give a formally ac¬ 
ceptable fit to the late-time radio emission, it is roughly consis¬ 
tent with the most dominant features in the data. Fig. 5 shows the 
Hubble-flow model light curves that provide the best fitting to the 
late-time radio data (where we have added eight data points pub¬ 
lished in Krauss et al. (2011) from days 141 and 175 to our mea¬ 
surements to help constrain the turnover of the high-frequency light 
curves). Taking the nova to have a distance of 6 kpc and using the 
expansion velocity measured by Yamanaka et al. (2010) on the day 
of discovery as the maximum velocity v^bs = 1500 km s“*, ejected 
at t = 0, and allowing other parameters to vary, our best fitting 
model has an ejecta mass of M = 1.7 x 10“"* Mq, a temperature 
of r = 1.1 X 10** K, and a ratio of minimum to maximum veloc¬ 
ity f =v,7vobs = 0.17 (see Fig. 5). Given that the ejecta are clearly 
not spherically symmetric (see Fig. 3), and that the evolution of the 
spectral index tttrans from 0.6 to 0.3 suggests a flattening of the den¬ 
sity profile over time, it is not too surprising that the Hubble-flow 
model does not provide a formally acceptable fit. Overall, however, 
the fitting of the late-time radio data with a Hubble-flow model in¬ 
dicates that if the ejecta were approximately spherical, then roughly 
2 X 10“"* Mq was ejected at a temperature of 10"^ K. 

Comparing model light curves for a range of ejecta masses, 
temperatures, and shell thicknesses provides a measure of the un¬ 
certainty on these parameters. Fig. 6 suggests that the best fitting 
parameters are accurate to within at least a factor of 2. In terms of 
the effect that modifying these parameters has on the model light 
curves, changing the distance would have influenced the maximum 
flux density at each frequency, to which it is inversely related. In¬ 
creasing the temperature of the shell increases the resulting flux 
density at all times, and causes the turnover from optically thick to 
thin to occur earlier. Increasing the maximum velocity also causes 
this transition to occur at earlier times. The ratio of the minimum 
and maximum velocity, f, relates to the flux density and the decline 
of the light curve after peak - a higher ratio corresponds to a ge¬ 
ometrically thinner remnant, which both raises the maximum flux 
density and causes a swifter decline in flux density after this peak. 


Whereas the late-time emission is consistent with originating 
in a thermal remnant with T ~ 10'‘ K and a strong resemblance 
to a Hubble flow, the early-time emission is not (see Fig. 7). In 
fact, the T ~ 10"* K ejecta were only revealed once the early-time 
flare faded. The brightness temperature at early times was strongly 
dependent on both frequency and time, with a maximum Tj of be¬ 
tween ~ 6x 10“* K at 36.5 GHz and ~ 4x 10® K at 1.4 GHz. Bright¬ 
ness temperature is an expression of surface brightness and rep¬ 
resents a lower limit on the physical temperature of any thermally- 
emitting material; for the flare to have been thermal, it would there¬ 
fore have had to have been emitted by plasma with T > 4 x 10® K. 

During the early-time flare, however, the X-ray flux was too 
low for the radio flare to have been due to hot (X-ray emitting) 
plasma. For the radio flare to have been entirely thermal, it would 
have had to have been produced by plasma that was optically thick 
at 1.9 GHz, had an electron temperature T a; 4 x 10® K, and had a 
size scale of 3 x lO^”* pc ( i 5 oo\ms~‘ )(4j^)’ could in prin¬ 

cipal also have been produced by free-free emission from hotter 
plasma that subtended a smaller solid angle on the sky. Krauss et al. 
(2011) detected X-ray emission from plasma with T > 2x10^ K on 
day 41 after the start of the eruption. For ejecta with T = 2 x lO’ K 
to be optically thick, the free-free optical depth would have to have 
been greater than one: 


Tff = 3.0 X 10 


^2x 10’ K 
This requirement corresponds to 
T 


1.9 GHz 


EM,,, >1. (4) 


f 


ntdl > 3.3 X 10' 


2 X 10’ K 


1.9 GHz 


pc cm. 


(5) 


A hot plasma with a temperature five times the maximum Tj, could 
in principal have produced the radio flare if it was optically thick 
at 1.9 GHz and subtended a solid angle that was 20% that of the 
entire ejecta. A uniform-density spherical shell of that size, with 
shell thickness I, would therefore require an electron density of 
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Figure 5. The late time observations of V1723 Aql between 2011 March and 2014 March, presented here for the first time, with eight previously published data 
points from days 2011 January 31 and 2011 March 5 used to constrain the high frequency light curves. Model indicates light curve of an expanding spherical 
thermal shell with a density profile of 1/r-, mass M=1.7 x 10““^ M©, ejecta temperature of T=l.l x 10"^ K, maximum velocity of Vobs=1500 km s,“^ distance of 
d= 6.0 kpc, and ratio of minimum to maximum velocity ^ = v//Vobs = 0.17. 
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( 6 ) 


to be optically thick. Such hot plasma would produce an X-ray 
emission measure, EM^ = f n^dV (where V is the emitting vol¬ 
ume), of approximately 8 x 10^* cm“^ (independent of £), and 
(unless the shell was extremely thin) have a mass of a few times 
10“'’ Mq. On day 41, however, Swift/XRT detected X-rays with 
an £Mx of just ~ 10^® cm“^, indicating that there was not enough 
2 X 10^ K plasma to produce the observed radio flare. Even gen¬ 
erating a model of the expected X-ray emission and assuming the 
highest reasonable absorbing column allowed by the X-ray spec¬ 
tral fit (3 X 10^^ cm““), the observed X-ray flux was still an order of 
magnitude lower than would be required to explain the radio flare 
with thermal emission. The natural alternative to thermal emission 
is radio synchrotron emission from electrons accelerated in a shock. 

The inconsistency between the radio flux density at 1.9 GHz 
and thermal radio emission expected from the 2 x 10’ K X-ray 
emitting plasma is independent of the distance to V1723 Aql. If 
we were to allow the temperature of the hot plasma to drop below 
2 X 10’ K (despite the constraints from Krauss et al. (2011)), then 
the evidence for synchrotron emission would in principle depend 
on the distance to V1723 Aql. However, for the distance to be small 


enough to explain the radio flare with thermal emission, the outflow 
velocities become inconsistent with constraints from optical spec¬ 
troscopy. Therefore, the combination of radio. X-ray, and optical 
observations appear to support a synchrotron emission mechanism. 


4.2.2 Failure of thermal shock emission model 

To further assess whether the early-time radio flare in 
V1723 Aql was thermal or non-thermal (synchrotron) in origin, we 
consider a one-dimensional, spherically-symmetric model for inter¬ 
nal shocks in novae, as presented in Metzger et al. (2014, hereafter 
M14). Consider a fast outflow of velocity Vf > 10’ km s“' that col¬ 
lides with a slower shell of velocity v^ < Vf. For simplicity, assume 
that the mass of the slow shell is comparable to or lower than that 
of the fast material and is emitted at time r = 0. The fast outflow is 
emitted a time At after the slow shell, resulting in a collision radius 


VsVfAt .3 Vs / At \ 

Rc ~ -6 X 10 Vf,8 — I-^ I 
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Av \ week ) 


and time 
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Figure 6. Variations of the Hubble Flow model compared to the observed flux densities. Unless otherwise stated, the parameters for the model show distance of 
6 kpc, maximum velocity of 1500 km s“^, temperature of 1.1 x lO"^ K, mass of 1.7x lO^'^M©, and ^ = 0.17. The leftmost column shows variation in temperature, 
with the top panel showing a temperature of 3 x 10‘^ K and the bottom showing a temperature of 5 x 10^ K. The middle column shows variation in mass: the top 
panel shows a model with mass 2.5 x lO”'^ A/©, and the bottom panel shows 1.2x lO”"^ M©. The right column shows valuation in shell thickness, varying the ratio 
of minimum to maximum velocity The top panel shows a ratio of 0.4, and the bottom shows a ratio of 0.08. 



Figure 8. Observed flux densities and model flux densities in radio for a 
thermal shock. Assuming a distance of 6 kpc, the model uses a fast wind 
with an initial velocity vy = 2100 km s“^ mass Mw = 0.8 x 10“‘^ M©, and 
a characteristic duration of 12 days. The wind collides with a slower moving 
circumbinary shell with initial velocity v^ = 1300 km s“^ and mass = 
1 X 10“‘* Mq. This results in a post shock velocity of Vobs =1500 km s“^ 
and total nova ejecta mass as = 1.8 x lO^'^M©, consistent with our late-time 
thermal model. 



Time Since Outburst (days) 


where Av = Vf - Vs and Vf,8 = Vf/10^ cm s~K 

The predicted radio flare following this collision would be 
powered by the forward shock that emerges from the slow shell 
at approximately the time of the collision, i.e. tc ~ 50 days. If 
both shells are ejected over a timescale less than the duration of 


the peak optical activity, i.e. At = (Av/vf)fc < 2 weeks, then 
the relative velocity between the shells must be relatively low, 
< Vf(A?/rpk) » 0.3vf, corresponding to a collision radius 
Vs?c < Vf?c 10‘^ cm for Vf 2 x 10^ km s ~^. 

The forward shock heats gas to a temperature (M14, eq. 31), 

1.4x lofl— —-] 1.3x 10'^(AL] K (9) 

\l03kms-V \0.3vfj 

corresponding to an X-ray energy of kT ~ 0.5 keV for characteris¬ 
tic values Av 0.3vf and Vf ~ 2 x 10^ km s“'. For purposes of an 
estimate, here we have assumed that the slow shell is of comparable 
or lower density than the fast outflow. 

For n sj 10^ cm“^ and kT x 0.5 keV the cooling time of the 
gas is approximately one week for solar abundances, and even less 
for Xz ~ 0.1 as assumed below. This suggests that the shock is 
likely to be radiative, in which case the resulting X-ray luminosity 
is given by 



where n is the (unshocked) density of the slow shell at the time 
of the collision. That the observed X-ray luminosity of Lx ~ 3 x 
10^^ erg s“' near the radio peak (Krauss et al. 2011) is several orders 
of magnitude smaller than the fiducial value in equation (10) would 
require strong absorption of these soft X-rays by neutral gas ahead 
of the shock. Indeed, the bound-free X-ray optical depth ahead of 
the shock (M14, their eq. 50) 
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Figure 7. Flux densities (top panel) and brightness temperatures (bottom panel) for V1723 Aql between 2010 September and 2014 March. The 
dashed lines indicate the expected thermal emission from an expanding spherical shell with a density profile proportional to 1/r^, mass M=1.7 x 
10“"^ Mq, ejecta temperature of T=l.l x lO'^ K, maximum velocity of Vobs=1500 km s“^ distance of d = 6.0 kpc, and ratio of minimum to maximum 
velocity ^ = v//Vobs = 0.17. The brightness temperature, which is an expression of surface brightness, is highly frequency dependent and peaks 
during the early-time flare at all frequencies. The great difference between the brightness temperatures during the early-time flare and the late-time 
peak highlights the physically distinct nature of the early and late radio brightenings. 
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exceeds unity for n > 10’ cm“’, where Xz = 0.1 is the assumed 
mass fraction of CNO elements. This density corresponding to a 
minimum mass of the slow shell of Msh > 4n(RlA)nmp ~ lO^'* Mq 
for Rc ~ 10*^ cm. For values of n helow which < 1, in equa¬ 
tion (10) still exceeds the observed X-ray luminosity by at least an 
order of magnitude; this may suggest that Vy g < 1, or that a frac¬ 
tion of the shock power instead emerges at UV/soft X-ray energies 
below the Swift bandpass. 

The maximum radio brightness temperature of radiative 
shocks at frequency v is given by (M14, eq. 65), 


3 X lOMn K ~ 10® K « > «th 

LTo.vJ r-lOS-lO^K 

To,v = 10® rij Vf,8 K » 4 X 10® vf (^) K n< «th, 
where Vm = v/10 GHz, iij = nj IQ’ cm,“® and 
/ A \~ ^ 

nth » 4 X 10’ ( 1 Vf;g v’o cm^® (13) 

defines the pre-shock density above which the layer of X-ray ion¬ 
ized gas ahead of the shock is optically thick to free-free absorption, 
where a temperature of 2 x 10"* K is assumed for the pre-shock gas. 
This value, appropriate to the photo-ionized layer slightly ahead of 
the shock (Metzger et al. 2014), may be slightly higher than the 
values expected for the bulk of the ejecta later, which our Hubble 
flow modeling finds is closer to 1 x 10** K. 

The observed peak brightness temperature of V1723 Aql (see 
Fig. 7) exceeds the range given by equation (13) by more than order 
of magnitude, strongly suggesting a non-thermal origin for most of 
the observed radio emission. Furthermore, equation (13) was de¬ 
rived under the assumption that free-free emission is the sole source 
of cooling behind the shock. If the effects of line cooling are in¬ 
cluded, Equation (13) greatly overestimates the brightness temper¬ 
ature, further enhancing the tension with thermal emission models. 

It is difficult to produce a thermal model that fits the early 
time radio data and is also consistent with the X-ray and late ra¬ 
dio constraint. In Fig. 8, we applied the M14 model to the early 
light curve, modeling the shock as the collision of a fast wind 
with = 2100 km s“' and mass 0.8 x 10“"* Mg and a shell with 
V., = 1300 km s“* with mass M^h = 1 x 10^“* Mg. We assumed 
that the shells were released around t = 0, with a characteristic du¬ 
ration of the wind of 12 days, motivated by the observed optical 
decline rate. These parameters result in a post-shock velocity of 
Vobs = 1500 km s“* and a total ejecta mass of M ~ 1.8 x 10“"^ Mg, 
consistent with our late-time thermal model. Taking n oc 1/r’ and 
distance of 6 kpc, we were unable to reproduce the high radio 
flux densities over the course of the shocks. We conclude that syn¬ 
chrotron emission must therefore have at least partially powered 
the early-time flare. 


4.2.3 Synchrotron emission from particles accelerated in shocks 

The failure of thermal emission to account for the early-time 
radio flare suggests that the flare was dominated by synchrotron 
emission from particles that were accelerated in shocks. For one 
thing, shocks are common in classical novae. Evidence for shocks 
is provided by X-ray emission from plasma with temperature above 
10’ K (e.g. Mukai et al. 2008). Such temperatures arise naturally 
behind shock with speeds on the order of 1000 km s“' (see Equa¬ 
tion 9). Furthermore, the recent identification of classical novae as 
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a new class of y-ray sources shows that not only are shocks ubiq¬ 
uitous in novae, but that relativistic particles are as well (Acker- 
mann et al. 2014; Metzger et al. 2015a). A natural acceleration 
mechanism for these particles is Fermi acceleration in the shocks. 
In fact, radio observations have illustrated the existence of radio 
synchrotron emission in a normal nova - in Fermi-detected nova 
V959 Mon, Chomiuk et al. (2014) identified the location of y-ray 
producing relativistic electrons as the region where a fast outflow 
collided with, and was shaped by, a dense equatorial torus. Al¬ 
though V959 Mon did not produce a strong early-time radio flare, 
its early-time radio emission was consistent with being synchrotron 
(Chomiuk et al. 2014). Whereas existing observations do not con¬ 
strain the inclination of V1723 Aql, V959 Mon is viewed very 
nearly edge on (Chomiuk et al. 2014) - it is thus possible that some 
early non-thermal emission may have been obscured by the opti¬ 
cally thick equatorial torus. Moreover, while the early-time flare in 
V1723 Aql suggests stronger non-thermal radio emission than that 
seen from V959 Mon, VLA observations of y-ray producing nova 
V1324 Sco (Nova Sco 2012) revealed a double peaked radio light- 
curve similar to that of V1723 Aql (Finzell et al., in preparation). 

Moreover, synchrotron emission is physically plausible for 
V1723 Aql. We expect the magnetic field behind the shock to 
be amplified to a characteristic value (e.g. Caprioli & Spitkovsky 
2014) 


B = 


{6neBmpnAwf’’^ 

0-054!i(] 


H0’( 



(14) 


where eg = 10 “’eB _2 is the fraction of the post-shock thermal en¬ 
ergy placed into turbulent magnetic fields. The characteristic syn¬ 
chrotron frequency from an electron of Lorentz factor y,, is given 
by 



For typical densities in the radio-producing nova shocks, n ~ 
10’ cm,“® the production of radio synchrotron thus requires rela¬ 
tivistic electrons with y^, corresponding to energies of a few hun¬ 
dred MeV. Such electrons may be directly accelerated at the shock 
(leptonic scenarios; see, e.g.. Park, Caprioli & Spitkovsky 2015; 
Kato 2015), or they may be created as electron/positron pairs pro¬ 
duced by the decay of charged pions (hadronic scenarios). The Lar- 
mor radii of the radio emitting particles (y^ of a few hundred) is 
extremely small compared to the size of the shock, so particles can 
easily remain within the acceleration region. As a comparison, su¬ 
pernovae remnants, which typically have post-shock densities on 
the order of 0.1 - 10 cm,“® and velocities of a few thousand km s“' 
(e.g., Williams et al. 2014; Metzger et al. 2015b), using equation 14 
would have magnetic field strengths scaling roughly as 50 pG, not 
atypical for post-shock remnants (as in, e.g., Ressler et al. 2014). 
More detailed modeling of the non-thermal radio emission in no¬ 
vae (Vlasov et al., in preparation) find acceptable fits to the peak 
brightness temperature for assumed microphysical shock parame¬ 
ters (Ce ~ 10“® - 10“’,eg ~ 10“’) as inferred from modeling super¬ 
nova remnants (e.g., Morlino & Caprioli 2012). 

The properties of the radio flare are also consistent with a syn¬ 
chrotron contribution. The expected time scale for the decay of 
any synchrotron contribution is set by the size scale of the rem¬ 
nant divided by the expansion speed (Taylor et al. 1989), which for 
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V1723 Aql is around 40 days, as observed. Krauss et al. (2011) did 
not detect linear polarization at 5.5 GHz on day 48, with a 3cr up¬ 
per limit of 39 /rJy (or 1%). Similarly, on day 57, we detected no 
linear or circular polarization averaging across 1 GHz of bandwidth 
centred at 28.2 GHz, with 3cr upper limits of 0.42 mjy (9%) for cir¬ 
cular polarization and 0.63 mJy (6%) for linear polarization. Given 
that the magnetic field geometry was likely to have been complex, 
the lack of strong polarization is not a concern for the synchrotron 
picture. That the radio spectrum rose with frequency during the ra¬ 
dio flare indicates that the synchrotron emission was either self- 
absorbed or suffered from free-free absorption. Since synchrotron 
self absorption would require the synchrotron-emitting region to 
have been extremely compact (i.e., have had a solid angle multiple 
orders of magnitude smaller than that of the remnant), whereas a 
reasonable amount of 10** K ionized gas could have led to observ¬ 
able absorption, we find free-free absorption more likely. 


4.3 Nova Shocks in the context of bipolar flows and y-rays 

The perpendicular structures evident in Fig. 3 are reminis¬ 
cent of the perpendicular flows thaf Chomiuk et al. (2014) found 
to be responsible for the generation of shocks and gamma-rays in 
nova V959 Mon. In V959 Mon, radio observations showed that 
strong internal shocks produced synchrotron emission coincident 
with the y-ray detection (Chomiuk et al. 2014). Radio images of the 
ejecta showed rapidly expanding bipolar flows, which faded over 
the course of sixteen months to reveal an inner structure that colli¬ 
mated the outer flow (see Chomiuk et al. 2014, Fig. 2d). Our image 
from day 1281 (see Fig. 3, panel b) reveals a similar inner struc¬ 
ture that may be a shell with an equatorial density enhancement. 
The evolution of the radio spectrum of V1723 Aql may also indi¬ 
cate that the inner portions of its ejecta had a flatter density profile 
than the outer regions, as one might expect if an inner structure was 
more confined fo the orbital plane (as in V959 Mon). Furthermore, 
the resolved image on day 832 showed two distinct components 
aligned roughly perpendicular to the axis of greatest elongation, 
which could be indicative of a limb-brightened collimating torus. 

By associating V1723 AqTs multiple, distinct radio structures 
with different velocities, we identified flows that likely collided 
and led to shocks and particle acceleration. While optical spec¬ 
tra revealed the presence of a fast flow with maximum radial ve¬ 
locity Vobs = 1500 km s“* (Yamanaka et al. 2010), there is opti¬ 
cal spectral evidence for slower moving ejecta as well. By using 
the break frequencies of the radio spectrum as the nova shell tran¬ 
sitioned from optically thick to optically thin, we estimated that 
the slowest ejecta had an expansion speed of V; between 200 and 
400 km s“*(d/6kpc)(ri,/10‘* K)“*^^. The two components that were 
visible on day 832 had a separation of 70 mas on the plane of the 
sky, which, using a distance of 6 kpc, implies that each component 
had a velocity of roughly 440 km s“' if they were ejected at t = 0 
- consistent with our estimate for the innermost material. Addi¬ 
tionally, Nagashima et al. (2013) modeled optical emission lines on 
days 7 and 11 with two overlapping gaussian components with cen¬ 
tral velocities of -190 km s“* and 360 km s“' relative to line of sight, 
also supporting a non-spherical distribution of ejecta with slow 
flows on the order of a few hundred km s“'. Finally, the expected 
temperature behind a shock between two flows with velocities of 

440 km s“* and 1500 km s“* is 1.4 x lO’ (tttttt— r) K (from 
equation 9), which is consistent with observational constraints on 
the temperature of the X-ray emitting plasma (Krauss et al. 2011). 
That we measure a density profile close to \jr^ from the transi¬ 


tioning ejecta indicates that the bipolar outflow likely had a wide 
opening angle. A model of V1723 Aql as a wide, biconical flow 
with collimating torus expanding at a few hundred km s“' is con¬ 
sistent with imaging, optical and radio spectroscopy, the detection 
of X-ray emission, and the evolution of the radio light curve. The 
presence of more complex substructures or additional shocks could 
account for other variations in the radio light curve. 

Fast biconical flows and slow equatorial tori may be more 
common than previously thought. For example, a change in the 
PA of maximum elongation like that in V959 Mon and V1723 Aql 
was also evident in Nova QU Vul, which had a double peaked light 
curve much like that of V1723 Aql (Taylor et al. 1987, 1988). Early 
imaging of QU Vul showed bipolar structure, later evolving into a 
more spherical shape - suggesting that as the outer bipolar emis¬ 
sion became optically thin, it revealed the central, more ring-like 
emission region (Taylor et al. 1987), as we saw in V1723 Aql. 

That V1723 Aql displayed evidence for both shocks and dust 
production (Nagashima et al. 2013) around the time of the radio 
flare raises the intriguing possibility that dust production in novae is 
related to shocks. Both X-rays and radio synchrotron emission indi¬ 
cate that shocks were energetically important in VI723 Aql roughly 
40±20 days into the eruption. Nagashima et al. (2013) argued based 
on optical observations that dust was produced around day 20. 
The finding by Metzger et al. (2014) that internal shocks gener¬ 
ate a cool, dense post-shock region suggests that the appearance of 
shocks and dust at approximately the same time in V1723 Aql may 
be more than a coincidence. 


5 CONCLUSIONS 

Radio and optical imaging of novae has shown evidence for 
the formation of jets, eg. RS Oph (Sokoloski et al. 2008); rings, eg. 
DQ Her (Slavin et al. 1995; Sokoloski et al. 2013); clumping, eg. T 
Pyx (Shara et al. 1997) and V458 Vul (Roy et al. 2012); collimating 
torus, eg. V959 Mon (Chomiuk et al. 2014); and other asymmetries. 
For V1723 Aql, the contrast between the early- and the late-time 
radio emission reveals a system of colliding flows in the ejecta. 
This finding has broader implications for the production of y-rays 
and the shaping of the ejecta from novae. 

(i) Late-time radio emission from V1723 Aql was roughly consistent 
with a freely expanding thermal plasma with a well-defined inner 
edge, a temperature of 10"* K, an ejecta mass of approximately 2 x 
10“"* Mo, a ratio of minimum to maximum velocity ^ =v,7vobs ~ 
0.2, a maximum velocity of 1500 km s“', and a distance of 6 kpc. 

(ii) The early-time radio flare was most likely due to a combination 
of thermal and non-thermal synchrotron emission emanating from 
the shocks that produced the observed X-ray emission. Thermal 
shocks alone fail to reproduce the high brightness temperature and 
weak X-ray emission. 

(iii) The shocks that generated the early-time radio flare did not cause 
either the early optical light curve or the late-time radio evolution to 
significantly diverge from the expected behaviour of a standard fast 
nova. Relatively few novae have such comprehensive early-time ra¬ 
dio observations, and thermal X-ray emission from shocks in novae 
is known to be common (Mukai et al. 2008) - therefore, detectable 
radio emission from early shocks could be more ubiquitous than 
previously realized. 

(iv) It has been established that y-ray production in the first few weeks 
of nova eruptions may be common (Ackermann et al. 2014; Met¬ 
zger et al. 2015a). y-rays are presumably the result of particle ac¬ 
celeration in shocks, early radio flares - as in V1723 Aql - may 
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provide an important diagnostic of the shocks and colliding flows 
that produce y-rays. 

(v) Although the morpho-kinematic structure of the ejecta from 
V1723 Aql is not tightly constrained, its behaviour is consistent 
with that of a biconical flow loosely collimated by a slow equato¬ 
rial torus, as in V959 Mon (Chomiuk et al. 2014). This hypothesis 
is supported by the existence of the perpendicular structures in the 
resolved radio images (see Fig. 3). Moreover, the shock speed de¬ 
rived from the difference between the fast radio flow (1500 km s“*) 
and the slow radio flow (knots moving at ~400 km s“') is consis¬ 
tent with that inferred from the temperature of the post-shock X- 
ray emitting plasma. The flattening of the density profile seen from 
the variation in tttrans could be a consequence of the material being 
slightly more concentrated towards the equatorial plane in the inner 
region of the remnant (see Fig. 4). 

Radio observations are a powerful diagnostic for the colliding 
flows and shocks that are crucial for y-ray production and shap¬ 
ing of the ejecta from novae. The radio emission from V1723 Aql 
shows that the collision of fast and slow flows resulted in a shock 
which, although leaving the late-time behaviour of the radio light 
curve mostly unaffected, produced both thermal and non-thermal 
emission. With the recent detections of classical novae in y-rays, 
the question of relativistic particles in novae becomes even more 
important. Chomiuk et al. (2014) found the morphology of the 
shock-producing structures in V959 Mon; our work here extends 
those findings to other classical novae. 
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